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The highly diluted antiferromagnet Mno.35Zno.65J2 has been investigated by neu- 
tron scattering in zero field. The Bragg peaks observed below the Neel temperature 
(Tjv ~ 10.9 K) indicate stable antiferromagnetic long-range ordering at low temper- 
ature. The critical behavior is governed by random-exchange Ising model critical 
exponents (y ps 0.69 and 7 « 1-31), as reported for Mn x Zni- x F2 with higher x 
and for the isostructural compound Fe x Zni- x F2. However, in addition to the Bragg 
peaks, unusual scattering behavior appears for \q\ > below a glassy temperature 
T g ~ 7.0 K. The glassy region T < T g corresponds to that of noticeable frequency 
dependence in earlier zero- field ac susceptibility measurements on this sample. These 
results indicate that long-range order coexists with short-range nonequilibrium clus- 
ters in this highly diluted magnet. 



Diluted uniaxial antiferromagnets have been extensively studied as physical realizations of theo- 
retical models of random magnetism Q] , including those pertaining to percolation phenomena [^],[| . 
For three dimensions (d = 3), two of the most extensively studied examples are the rutile com- 
pounds yLn x Zn\—xF% and Fe x Zri\- x Fi. These two systems differ effectively only in the strength 
and nature of the anisotropy, providing a unique opportunity to explore the role of anisotropy in 
the ordering of dilute magnets at low temperature. In Mn x Zn\- x Fi the anisotropy is dipolar in 
origin Q. In Fe x Zn\- x Fi the anisotropy is an order of magnitude greater for x = 1 because of 
the additional crystal field contribution || . In many experiments with the magnetic concentration, 
x, well above the percolation threshold concentration x p = 0.245 J2j, the behaviors for H = are 
qualitatively similar for Mn x Zn\- x Fi and Fe x Zn\- x F2- Antiferromagnetic (AF) long-range order 
(LRO) at low temperatures and characteristic random-exchange Ising critical behavior have been 
observed in the Fe x Zn\^. x Fi compounds for x > 0.31 ||. Similar random-exchange Ising model 
(REIM) behavior is found in the Mn x Zn 1 _ x F 2 system for x > 0.4 0. 

For small fields applied parallel to the uniaxial direction and reasonably small magnetic dilution, 
the diluted antiferromagnet in a field (DAFF) is expected to show critical behavior belonging to 
the same universality class as the random-field Ising model (RFIM) for the ferromagnet, the latter 
being the model most used in simulations Jl0[ . Indeed, for all measured samples of both systems 
for which the REIM character is found at H = 0, the application of a small field parallel to the 
easy axis generates critical behavior compatible with the predicted Q REIM to RFIM crossover 
scaling. In spite of the evidence supporting the DAFF as a realization of the RFIM, some non- 
equilibrium features inherent to DAFF compounds and also the newly explored field limits (Ii]-[l| 
of the weak RFIM problem in d — 3 make the nature of the phase transition at T C (H) still a matter 
of considerable controversy [|l4|-[l6) . 

Under strong random fields (corresponding to large H) and also close to the percolation threshold, 
the phase diagrams of DAFF's have proven to be much more complicated than originally anticipated. 
For large H, AF LRO is predicted to become unstable ^J§. The generation of strong random fields 
induces |l8| , |l9| a glassy phase in the upper part of the (H,T) phase diagram of d = 3 Ising DAFF's. 
The equilibrium boundary, 

T eq {H) = T N -bH 2 -C eq H 2 't (1) 

above which hysteresis is not observed, has a convex shape at high H (</> > 2), instead of the concave 
(<f) = 1.4) curvature seen at low field (where REIM to RFIM crossover occurs). This change of 
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curvature in T eq vs. H was first observed pC| by magnetization measurements in Feo.31Zno.a9F2- 
Faraday rotation H and neutron scattering J2l[ | experiments on a sample with the same x, confirmed 
the REIM to RFIM crossover scaling at low H and the lack of stability of the AF LRO at large 
H, giving way to a random-field induced glassy phase in this highly diluted compound. Recent 
magnetization measurements indicated that similar structure in the phase diagram exists at very 
high fields for samples with higher values of x |Tl] . |l2]]. At a still higher concentration the low 
temperature hysteresis observed for x < 0.8 is absent Tp2]. 

The magnetic features observed at large H in samples of Fe x Zn\^ x F2 in the concentration range 
0.3 < x < 0.6 contrast with the behavior in weakly anisotropic system Mn x Zn\- x Fi for intermediate 
x, where a strong H induces a spin- flop phase [p3| . This distinct behavior may be solely a consequence 
of the stronger Ising character of the former system. In the strong dilution regime (x « x p ), a number 
of magnetic features lead us to distinguish between these two systems, as well. For x < 0.27, no 
long range order |m| is observed in Fe x Zn\- x F2- Typical spin-glass behavior was found [^5| in 
a sample with x — 0.25, although recent works f26|-|28|| suggest non-critical dynamics for x close 
to x p in this system. Close to the percolation threshold even a minute exchange frustration is a 
suitable mechanism [^9| for the spin-glass phase in Fe x Zn\~ x F2, as supported by local mean-field 
simulations |30|. For Ising systems, it is also expected that the dynamics even at zero field should 
be extremely slow ]3l| . In Mn x Zni- x F2, ac susceptibility measurements indicates a spin-glass 
clustering at low temperatures for samples with Mn concentrations 0.2 < x < 0.35 |3^J&|. Earlier 
neutron scattering studies |54| suggest, however, that at H = the termination of the line of the 
AF-paramagnetic (P) continuous phase transition occurs at T = at x — x p in stark contrast to the 
behavior of Feo.2bZno. 75-^2- In light of this contrast, the influence of the frozen spin-glass clusters 
on the stability of the AF LRO for x close to x p in Mn x Zn\^ x F2 is an important question that 
motivated the present work. The dipolar anisotropy of this weakly anisotropic system is expected to 
become random in strength and direction as x decreases, in contrast to the x-independent single-ion 
anisotropy of Fe x Zn\- x F2- In the case of Mn x Zn\- x F2 under strong dilution, the application of 
the results from numerical simulations fl30] ] applied to Ising systems is of course not warranted. Any 
differences observed in this system and the Feo.31Zno.69F2 must certainly reflect the difference in 
anisotropy and this may give a window to the understanding of the general phase diagrams for dilute 
anisotropic antiferromagnets in applied fields. 

In this study we performed zero-field neutron scattering experiments in Mno.35Zno.05F2 to verify 
the existence of a stable long-range ordered antiferromagnetic phase below a critical temperature 
Tjv ~ 10.9 K, where REIM critical exponents v « 0.69 and 7 k, 1.31 govern the behavior and to 
investigate the dynamic features of the system at low temperature. An unusual scattering behavior 
appears, for \q\ > 0, below T g 7.0 K, corresponding to the region where earlier ac susceptibility 
studies |33| indicated a noticeable frequency dependence in the real part of the susceptibility, in the 
absence of external field. The results indicate that long-range order coexists with non-equilibrium 
clusters in this highly diluted system. 

The neutron scattering experiments were performed at the Oak Ridge National Laboratory using 
the HB2 spectrometer in a two-axis configuration at the High Flux Isotope Reactor. We used the 
(002) reflection of pyrolitic graphite to monochromate the beam at 14.7 meV. The collimation was 
60 minutes of arc before the monochromator, 40 between the monochromator and sample, and 40 
after the sample. A pyrolitic graphite filter reduced higher-energy neutron contamination. The 
c-axis of the crystal was vertical and parallel to the applied field. A small mosaic was observed 
from the Bragg peak scans at low temperature, with roughly a half-width of 0.2 degrees of arc or 
0.0035 reciprocal lattice units (rlu). The mosaic was incorporated into the resolution correction 
by numerically convoluting the measured resolution functions, including the mosaic, with the line 
shapes used in the data fits Most of the scans taken were (1 q 0) transverse scans. For simplicity, 
the line shapes used in the fits to the data are of the mean-field form 

S{q) = ^^ + M 2 s 5{q) , (2) 

where k = l/£ is the inverse fluctuation correlation length and M s is the Bragg scattering from the 
long-range staggered magnetization. The critical power-law behaviors are expected to be k = \t\ u , 
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X = Ak~ 2 = x^|£|~ 7 and M s — Molt] 13 , where M Q is non-zero only for t < 0. The exponents u, 
7, and j3 and amplitude ratios k + /k ~ and Xo + /Xo~ are universal parameters characterizing the 
random-exchange Ising model. 

The sample was wrapped in aluminum foil and mounted on an aluminum cold finger. A calibrated 
carbon resistor was used to measure the temperature. 

Transverse scans, taken after quenching to low temperatures (T = 5K) and subsequently heating 
the sample, are shown in Fig. 1. For clarity, the data for the range \q\ < 0.008 rlu, which spans the 
Bragg scattering component, are not shown. For the most part, the scans are quite consistent with 
what is expected for a phase transition occurring near T = \\K. However, a most unusual feature 
of the line shapes is evident in the data at the lowest temperature, T — 5 K. The broad line shape 
indicates a great deal of short range order present upon quenching. The short-range order is evident 
for the scans with T < 7 K. The scan at T — 6 K shows striking asymmetry as shown in Fig. 2. Since 
the scans were taken with increasing q from —0.19 to 0.19 rlu and each measurement took about 35 
seconds, the asymmetry is an indication that the short-range order is rapidly decreasing with time, 
i.e. the system is equilibrating. The slow relaxation for T < 7 K corresponds very well to the large 
frequency dependence observed using ac susceptibility in the same sample ||| for T < 7 K. 

A transition to antiferromagnetic long-range order is indicated by the presence of a resolution 
limited Bragg scattering peak which decreases sharply as T approaches T/v ~ H K. As T — > T/v 
from above, the width of the non-Bragg scattering component decreases and the q = intensity 
increases. Similarly, as T — > TV from below, the width decreases and the q = intensity increases. 
Such behavior is typical of an antiferromagnetic phase transition. To fit the data, we used the 
Lorentzian term in Eq. 1, convoluted with the instrumental resolution. The data for \q\ < 0.008 rlu 
were eliminated from the fits to the Lorentzian term to avoid Bragg scattering. The results of the fits 
yield re(T) and the staggered susceptibility x(T) = A/n 2 . The results for re are shown in Fig. 3, along 
with the expected random-exchange critical behavior |3(||37j as indicated by the solid curves with 
v = 0.69 and re+ / k~ = 0.69. The overall amplitude of the solid curves is adjusted to approximately 
follow the data. A clear minimum re ~ 0.017 rlu is observed in the fitted values near T/v, indicating 
significant rounding due to a concentration gradient in the crystal p8| ]. The gradient rounding is 
most likely the cause of the deviations of the data from the fit away from the minimum as well. 
Nevertheless, the present data are plausibly consistent with random-exchange critical behavior when 
the significant rounding due to the concentration gradient is taken into account. 

Results for the logarithm of x vs - T are shown in Fig. 4. The random-exchange behavior p6| , ^7| , 
with 7 = 1.31 and Xo + /Xo~ — 2.8 and with the overall amplitude adjusted to approximately fit the 
data, is shown as the solid curves. The maximum in the data and the systematic deviations from the 
fit are indications of a significant gradient in the concentration, as we discussed with respect to Fig. 
3. Again the data are fairly consistent with a concentration-rounded random- exchange transition to 
antiferromagnetic long-range order. 

The Bragg intensity, obtained by subtracting the fitted Lorentzian scattering intensity from the 
total q = scattering intensity, is shown vs. T in Fig. 5. Once again, the data are fairly consistent 
with a random-exchange |39| transition (fi = 0.35) near T = \\K represented by the solid curve. 
The nonzero Bragg component above T = 11K is probably attributable mainly to concentration 
gradient effects. The precise shape of the Bragg scattering intensity vs. T in Fig. 5 must not be 
taken too seriously, particularly at low T, since it is known that severe extinction effects distort 
the behavior by saturating the measured value Q. In addition, for T < 7 K, the sample shows 
nonequilibrium effects since it was quenched, as described above, and the magnitude of the Bragg 
scattering component might well be smaller than if the sample were in equilibrium. The large Bragg 
scattering component well below T/v, along with the minimum in re and maximum in x near T/v 
strongly indicate an antiferromagnetically ordered phase. 

Previous magnetization studies indicate a de Almeida-Thouless-like (AT) curve in the 

H — T phase diagram. The H — endpoint of this boundary coincides reasonably well with the 
antiferromagnetic phase transition observed with neutron scattering. 

In conclusion, we have shown neutron scattering evidence that this system, Mno,asZno. 65-^2; or- 
ders near T = 11 K in a way consistent with the REIM model. In addition, significant relaxation 
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takes place for T < 7 K. This is consistent with previous magnetization measurements and demon- 
strates that only part of the system orders with long range order when the system is quenched to 
low temperatures. This behavior is consistent with clusters coexisting with long range order below 
Tm- A similar glassy low temperature region has been identified f|pO 
Fcq. 31 Zno. 69-F2 using magnetization and dynamic susceptibility measurements. However, the broad 
line shapes that indicate the glassy behavior were not observed in Feo. 31^^0.69^2 with neutron 
scattering techniques |2l| . It is interesting that neutron scattering measurements at the percolation 
threshhold in Mno.25Zno.75F2 did not indicate any glassy behavior in contrast to Feo.25Zno.75F2- 
This should be investigated further. 
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FIG. 1. The logarithm of the neutron scattering intensity vs. q just above and below TV in zero field 
obtained after quenching the sample to T = 5 K and heating. 
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FIG. 2. The logarithm of the neutron scattering intensity vs. q for T = 6.00 K in zero field obtained 
after quenching the sample to T = 5 K and heating. The data were taken in sequence in increasing q with 
approximately 35 seconds per point. The asymmetry indicates that on this time scale the system is relaxing 
toward the behavior seen at higher T. 
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FIG. 3. k vs. T near Tjv. The solid curves represent the expected random-exchange critical behavior. 
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FIG. 4. The logarithm of \ vs - T near Tn- The solid curves represent the expected random-exchange 
critical behavior. 
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FIG. 5. The Bragg scattering intensity. The solid curves represent the expected random-exchange critical 
behavior, vs. T. 
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